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model,  analytical  expressions  were  obtained  for  responses  of  the  beam  lying 
over  a  ground  surface  and  subjected  to  vertical  loads.  These  expressions  were 
then  coded  into  the  computer  program,  2PARM.  The  results  were  compared  with 
computed  results  obtained  by  a  finite  element  method. 

The  study  has  maintained  that  this  model  and  computer  program  can  well 
predict  the  responses  of  the  beam  bearing  on  the  surface  of  an  elastic 
homogeneous  medium. 


PREFACE 


This  report  presents  investigations  and  results  of  utilizing  a  two- 
parameter  layer  soil  model  for  solving  soil-structure  interaction  problems. 

The  two-parameter  models  are  capable  of  coupling  the  responses  of  the  soil  to 
the  structure,  yet  are  simple  enough  to  b^  developed  as  analytical  expressions 
in  an  easy-to-use  computer  program. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


_ Multiply _ 

feet 

foot-pounds  (force) 

kips  (force) 

kip  (force)-feet 

kips  (force)  per  square  foot 

kips  (force)  per  square  inch 

pounds  (force) 

pounds  (force)  per  square  foot 
pounds  (force)  per  square  inch 


By 

0.3048 

1.355818 

4.448222 

1355.818 

47.88026 

6.894757 

4.448222 

47.88026 

6.894757 


_ To  Obtain 

metres 
metre-newtons  or  joules 
kilonewtons 
newton-metres 
kilopascals 
megapascals 
newtons 
pascals 
kilopascals 


SOIL-BEAM  INTERACTION  ANALYSIS  WITH  A  TWO-PARAMETER 
LAYER  SOIL  MEDIUM:  HOMOGENEOUS  MEDIUM 


PART  I:  INTRODUCTION 

One-Parameter  Versus  Two-Parameter  Models 

1.  In  the  analysis  of  a  beam  bearing  on  a  ground  surface,  the  soil 
medium  is  often  idealized  by  a  Winkler  one-parameter  model,  significantly 
simplifying  the  representation  of  the  ground  behavior.  This  idealization 
makes  it  possible  to  develop  analytical  expressions  for  the  beam-soil  system 
for  a  variety  of  loading  and  boundary  conditions.  Unfortunately,  the  soil 
response  at  one  point  is  uncoupled  from  the  displacement  at  other  points  in 
the  Winkler  model,  and  the  model  fails  to  reproduce  the  characteristics  of  a 
continuous  medium.  To  resolve  this  problem,  two-parameter  models  have  been 
developed  (Filonenko-Borodich  1940,  Pasternak  1954).  These  models  are  capable 
of  coupling  the  response  within  the  soil,  yet  simple  enough  to  permit  develop¬ 
ment  of  analytical  expressions  for  the  beam-soil  system.  This  is  accomplished 
with  more  or  less  the  same  degree  of  difficulty  as  that  associated  with  the 
Winkler  model. 

Parameter  Determination 

2.  The  parameters  of  the  model  can  be  determined  from  known  displace¬ 
ments  of  the  ground  surface  (Kerr  1976)  or  by  using  a  variational  method 
together  with  an  assumed  displacement  profile  of  the  ground  with  depth  (Vlasov 
and  Leonteve  1966,  Harr  et  al.  1969).  However,  since  the  ground  surface  dis¬ 
placement  and  the  displacement  profile  are  not  known  before  the  analysis,  it 
is  difficult  to  determine  the  parameters  for  the  analysis.  When  the  param¬ 
eters  are  defined  from  the  ground  displacements  for  a  particular  case,  they 
may  not  be  applicable  to  other  cases.  Nonhomogeneity  of  the  soil  medium  makes 
the  above  stated  difficulty  even  more  pronounced.  Thus,  in  the  present  form 
of  the  two-parameter  models,  it  is  not  only  inconvenient  to  determine  the 
model  parameters,  but  also  is  very  difficult  when  the  soil  profile  is  the  only 
information  available. 
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3.  The  objective  of  the  study  is  to  develop  an  economical,  analytical 
tool  for  the  response  analysis  of  a  beam  resting  on  ground  surface.  To 
achieve  this  objective,  a  simple  mechanical  soil  model  is  developed  and  used 
in  formulating  the  responses  of  the  beam  resting  on  a  ground  surface.  The 
developed  model  is  classified  as  a  two-parameter  model.  Contrary  to  conven¬ 
tional  two-parameter  models,  however,  the  parameters  of  the  new  model  can  be 
defined  from  the  information  of  the  soil  profile  alone  for  any  complex 
profile . 


PART  II:  TWO-PARAMETER  MODEL  FOR  LAYERED  MEDIUM 


4.  A  horizontal  soil  deposit,  underlain  by  a  rigid  base,  is  idealized 
as  a  system  made  of  elastic  vertical  columns  as  shown  in  Figure  1.  The  system 
is  assumed  to  be  subjected  to  vertical  loads  at  the  surface  and  to  be  under 
plane-strain  conditions  (e.g.,  the  columns  and  loads  are  infinitely  long  in 
the  direction  perpendicular  to  the  cross  section  shown  in  Figure  1).  The 
forces  induced  in  the  system  are  axial  forces  in  the  column  and  friction 
forces  along  the  vertical  sides  of  the  column.  For  an  infinitely  small  width 
of  the  beam  element,  dx  ,  the  axial  normal  stress  and  surface  shear  stress 
are  expressed  by,  respectively, 


o(x,z) 


3w(x , z) 

E  “ 3i 


T (x , z )  =  G' 


3w(x,z) 

3z 


(1) 


where 


a  and  t  =  shear  stress  along  the  sides  of  the  column  and  axial 
normal  stress  in  the  column 

w  =  vertical  displacement  of  the  column 

E'  and  G'  =  constants  of  the  column  associated  with  axial  force  and 
friction  force 

dx  =  cross-section  area  of  the  column 


P 


Figure  1.  Soil  medium  idealized 
by  a  system  of  elastic  columns 


5.  The  displacement  w(x,z)  in  the  system  is  generally  expressed  by 


w(x,z)  =  I  <J>  (z)W  (x)  (2 

n  n  n 

where 

4>  (z)  =  assumed  function  of  z 
n 

W  (x)  =  unknown  vertical  displacement 
n 

In  order  to  define  <J>n(z)  »  two  approaches  were  considered  in  the  study.  In 
the  first  approach,  orthogonal  functions  are  selected  for  ^(z)  which  can 
satisfy  the  homogeneous  boundary  conditions  at  both  ends  of  the  column.  The 
displacement  is  the  resultant  of  superpositions  of  various  shapes,  4>^(z)  » 
weighted  by  W^(x)  in  this  approach. 

Second  Approach  Preferred 

6.  In  the  second  approach,  a  piece-wise  linear  displacement  is  assumed 
along  the  column.  Thus,  the  column  is  divided  into  a  number  of  segments  and 
one  linear  function  is  assumed  within  a  segment.  The  second  approach  is  more 
convenient  and  computationally  more  efficient  than  the  first  approach.  The 
expressions  and  computer  program  presented  in  this  report  are  based  on  the 
second  approach,  while  the  expressions  based  on  the  first  approach  are  given 
in  a  separate  volume. 

7.  Inside  the  model,  a  horizontal  layer  with  thickness,  i  ,  is 
considered  as  shown  in  Figure  2.  The  column  element  in  this  layer  is  sub¬ 
jected  to  vertical  forces  as  also  shown  in  Figure  2.  The  work  done  by  the 
internal  and  external  forces  for  this  element  are  expressed  by,  respectively, 


wi --- 1 


aw(x,z)  J__ 

a  — 5i —  dx 


W  =  {p  (x)  w(x,0)  +  p  (x)  w(x,«,)}  dx  -  f  t  w(x,z)  dz 


e  ‘  u 


+  /"  i  +  ~  ax  w(x,z) 

J  dX 


where 

W.  and  W  =  work  done  by  internal  and  external  forces 
i  e 

p^  and  p^  =  vertical  forces  acting  on  a  unit  area  of  the  upper  and 
lower  surfaces  of  the  layer 

o  and  t  =  stresses  explained  in  equation  1 

Substituting  equation  3  into  the  equilibrium  conditions  expressed  by 

Wi  +  We  =  0  ’  t*ie  f°H°win8  equation  is  obtained: 


/  fx  w(x’z^  dz  dx  +  J  a  w(x,z)  dz  dx 


=  Pu(x)  w(x, 0)  dx  +  p ^ (x)  w(x,«.)  dx 


u_  dx 

1 

r 

t 

U  A 

r  (x)dx 


a.  Horizontal  layer 


P.  ( x  )dx 
b 


b.  Column  element  subjected 
to  forces 


Figure  2.  Horizontal  layer  and  column  element  In  the  layer 

subjected  to  forces 


T 1  .VVW^  \,K  gg  A  ■■  y 


8.  When  a  continuous  linear  displacement  is  assumed  within  the  layer 
considered,  w(x,z)  in  equation  2  can  be  expressed  by 


w(x, z) 


(<P  ( z )) 

u 


Vz) 


T  / 


Wu(x) 


W£(x) 


(5) 


where  W  and  VG  =  vertical  displacements  of  the  layer  at  upper  and  lower 
ends,  respectively;  and 


IVZ> 


(6) 


Substituting  equation  1  into  equation  4  and  using  the  relationship  given  in 
equation  5,  equation  4  is  rewritten  by 


KM 

T 

'Wu(x)' 

|pu(x) 

'  +  [kl]  “ 

W*(x) 

|p£(x) 

(7) 


where  W"  =  d^W/dx^  ;  [NL]  and  [K^] 
of  the  column  and  are,  resnect ively , 


are  defined  by  the  material  constants 


[kl 


T 

Id*  (z)/dz|  (d^  (z)/dz) 
E'(z)  U  U  dz 

(d^^(z)/dz )  ( d1^  (z ) /dz  ) 


(8) 


Assuming  linear  variation  of  G'  and  E*  with  depth  z  within  a  layer,  the 
above  integrals  are  evaluated  as 


■-\a 

‘.■'•'-I 


m 


9 


(9) 


G1  (OH 
12 


E'(O) 


i  f3  n  g'(Q£  n 

'Li  iJ  12  Li 

r  i  -n  +  e '  ( £)  r  i 

L-i  iJ  £  L-i 


where 


G'(O)  and  E'(O)  =  material  constants  of  the  columns  at  the  upper 

surface  of  the  layer  considered 

G ’  ( Jl)  and  E *  ( fi.)  =  material  constants  of  the  columns  at  the  lower 

surface  of  the  layer  considered 

9.  When  the  system  is  divided  into  N  horizontal  layers,  the 
equilibrium  and  compatibility  conditions  at  the  interfaces  between  the  layers 


lead  to 


where 


[N]  and  [K] 


- [N] { W" (x) }  +  [K] {W(x) }  =  (P(x)} 


N  x  N  matrices  constructed  from  [NL]  and  [KL]  as 


shown  in  Figure  3 


[W]  =  vector  containing  vertical  displacements  at  N 

interfaces,  in  which  the  first  location  is  assigned  for 
the  displacement  at  the  ground  surface 

[P]  =  vectors  containing  external  forces  applied  at  the  N 

interfaces,  in  which  the  first  location  is  assigned  for 
the  force  applied  at  the  ground  surface.  The  expression 
is  identical  to  the  governing  equation  of  commonly  used 
two-parameter  models  except  in  a  matrix  form 

10.  When  a  beam  is  lying  horizontally  over  the  ground  surface  and  is 
subject  to  vertical  forces,  the  behavior  of  the  horizontal  beam  is  governed  by 
the  equation  expressed  by 


El  W,  (x)  =  P (x)  -  P  (x) 
b  s 


where 


El  =  flexural  rigidity  of  the  beam 

W,  =  vertical  displacement  of  horizontal  beam 
b 

„IV  ,4  ..  4 

W,  =  d  W,/dx 

D  D 

P  =  soil  reaction  force  at  the  soil-beam  interface 
s 


I'/’ \  • 

•  ■  —  ■  —  » m  ,t»  .m.m  m.. ».  .».#  . 


Figure  3.  Construction  of  [N]  and  [K] 
from  [NL]  and  [K^] 

Using  equations  10  and  11,  equilibrium  and  compatibility  conditions  at  the 
soil-beam  interface  lead  to 

[El]  (WIV(x)}  -  [N]  {W" (x)  }  +  [K]  { W (x) }  =  { P (x) } 

where  [El]  =  diagonal  matrix  containing  El  at  the  location  (1,  1)  and 

zeros  at  other  locations;  and  { P }  =  force  vector  containing  the  applied 
vertical  load,  at  the  first  location  in  the  vector  and  zeros  at  other 
locations . 


PART  III:  BOUNDARY-VALUE  PROBLEMS 


11.  Idealizing  the  soil  medium  as  a  mechanical  model  as  presented  in 
Part  II,  the  analytical  expressions  for  the  responses  of  the  ground  and  the 
responses  of  the  beam  bearing  on  a  ground  surface  are  obtained  solving 
equations  10  and  12,  respectively.  The  loading  conditions  considered  are 
vertical  loads  applied  on  ground  surface  and  vertical  and  moment  forces 
applied  to  the  beam. 

Responses  of  Ground  Without  Beam 

12.  The  expression  of  the  response  for  any  distributed  load  can  be 
obtained  by  integrating  the  expression,  for  a  line  load.  The  solution  for  a 
line  load  can  be  obtained  by  solving  equation  10  with  { P }  =  {0;  and  treating 
a  line  load  as  one  of  the  boundary  conditions. 

General  solution 

for  homogeneous  equations 

13.  The  solution  for  equation  10  with  { P }  =  {0}  is  expressed  by 


(W(x))  -  Y  U„>  ♦  C2n  .-X"»  <dn) 

n=l 


=  [ A(+x) ] { C 1 }  +  [A(-x) ] { C  2 } 


where 


{Cj }  and  {C2} 
[A(+x) ]  and  [A(-x)] 


X  and  {d  } 
n  n 


=  vectors  containing  the  unknown  constants 

=  N  x  N  matrices  in  which  the  columns  are 

Anx,,  ,  ,  Anx,,  , 

e  { d  }  and  e  {d  } 
n  n 

=  n^  eigenvalue  and  eigenvectors  obtained  from 
the  following  characteristic  equation  obtained 
from  equation  10 


[ -A  [N]  +  [K] I  (d)  =  {0} 


12 


With  equation  13,  the  vertical  shear  force  along  the  vertical  plane  are 
obtained  as 


{Tg(x)}  =  [N] {W1 (x) } 

=  [N]{[A'(+x)]{C1}  +  [A'(-x)]{C2}} 


(15) 


where 


[A' (ix) ] 


£b  A(±x)] ;  and 


{T  (x) }  =  vertical  shear  force 
s 


Response  due  to  a  line  load 
applied  on  a  ground  surface 

14.  The  soil  medium  is  divided  into  left-hand  and  right-hand  sides  with 
respect  to  the  location  of  a  line  load  applied  as  shown  in  Figure  4.  Since 


p 


Left  |  Right 

I 

I 


7V/  y  /  /  ///// 

Figure  4.  Line  load  applied  on  the  surface 

the  soil  responses  should  diminish  with  distance,  the  displacements  and  forces 
induced  in  the  soil  medium  (Equations  13  and  15)  are  deduced  to  the  following 
expressions : 


13 


lA  m.  .1  m 


v  t  v  "  .  <  ~  v  \  -  v  w  v  -  w  - 


*  .*  *.«  wj*  _w 


Left-hand  side 


{W(x))L  =  [A(+x)]{C}l  and  {Tg(x)}L  =  t N  ^  t A  *  C+x )  ]  { C  >L 


Right-hand  side 


{W(x) }R  =  [A(-x)](C}R  and  {Tg(x)}R  =  [ N ] [ A ' y (-x) ] {C }R 

where  the  subscripts  L  and  R  indicate  the  left-nand  and  right-hand  sides, 
respectively.  The  boundary  conditions  at  x  =  0  require 


{W(0)L  =  {w(0)  }R 


T  (0)  }  -  {T  (0) }  =  {P} 

S  L  S  K 


15.  Substituting  equation  16  into  equation  17  to  determine  the  unknown 
constants,  the  displacements  of  the  soil  medium  are  expressed  by 


(W(x)  }  =  ^  Ca(- |x |)J  {y} 


where  {y}  =  vector  containing  the  numbers  located  at  the  first  column  of  the 
matrix,  [ [N] [A* (—0) ] ]  1  .  Thus,  the  displacement  of  the  ground  surface 
subjected  to  a  line  load  at  x  =  i  is  expressed  from  equation  18  as 


V>  -*p 


X  e"'"lx' 


where 


W  =  displacement  of  the  ground  surface  (e.g.,  the  value  located  at  the 
first  location  {W} 

Y  =  the  value  located  at  the  n^  location  in  {y} 


Response  due  to  a  trapezoidal 
load  applied  on  ground  surface 


16.  Using  the  expression  in  equation  19,  the  surface  displacement  due 
to  a  distributed  load  applied  at  the  surface  is  obtained  as 


Wj (x)  = 


(2 


If  the  distributed  load  is  a  trapezoidal  load  expressed  by  P(£)  =  a£  +  b  , 
equation  20  yields  the  following  expression: 


Wj  (x) 


N 


I 

n=  1 


A 

n 


+  B 

n 


+  D 

n 


(2 


where 


C 


0  and  D 


0 


for  x  £  x„ 


General  solution  of 
homogeneous  equations 

17.  In  order  to  obtain  the  responses  of  the  beam  bearing  on  the  ground 
surface,  equation  11  with  {P}  =  {0}  must  be  solved.  The  characteristic 
equation  obtained  from  this  homogeneous  equation  is 

[A4[EI]  -  A2  [n]  -I-  [K]  ]  {d}  =  {0}  (23) 

N  sets  of  four  complex  eigenvalues  and  two  complex  eigenvectors  are  found  i 

t  ll 

from  equation  23  in  general.  The  n  four  eigenvalues  are  +A^  ,  -A  ,  +A^ 

and  -A'  in  which  A  and  A'  are  conjugate  pair.  The  eigenvectors 
n  n  n 

corresponding  to  ±A^  and  ±A^  are  also  conjugate  pair  and  denoted  as  {d^ } 

and  fd'}  ,  respectively.  Thus,  the  solution  of  the  homogeneous  equation  | 

above  considered  is  expressed  by 


A 

<"<*>>  ‘  J  Cl„  e  "  '“n1  +  C2n 

n=  1 


-A  x 


A '  x 
n 


n  (d  }  +  C,  e  (d '  ) 
n  3n 


(24) 


-  A '  x 

+  C.  e  n  id’  } 
4n  n 


where 

first 


and 


,  C_  ,  C.  ,  and  C.  =  unknown  constants 
2n  in  4n 

fourth  terms  and  between  the  second  and  third 


Combing  between  the 
terms  in  equation  24, 


this  equation  can  be  rewritten  as 


(W(x)}  =  2  Cln  e  {cos  (8nx)  ^dn)  -  sin  (B^)  jdJJ 
n=  1 


+  C2n  e  n  {cos  (Snx)  {d;}  +  sin  (B^)  {d^}} 


(25 


-a  x 


l  A  / 

+  C  e  n  {cos  (S  x) 
3n  f  n 


{dn}+  sin  (6nx){d;JJ 


-a  x 


+  C.  e 
4n 


n  {cos  (Snx)  {d;|  -  sin  (Bnx)  {d^| 


=  [A(+x)  ]  {C}  }  +  [  B  ( +x )  ]  t'C  }  +  [  A  (-x)  ]  ^C3  }  +  [B(-x)]  {C^} 


where  (d^  1  =  {d^|  +  i  {d*j  !  the  nth  columns  in  [A(ix)]  and  [B(±x)J  = 

e  n  {cos  (rnx)  {d^|  i  sin  (8^)  {dj  and  e  n  cos  (6^)  |d*| 

i  sin  (8  x)  Jd>  ,  respec C ive I v ;  and  A  =  a  +  iB 
n  J  n  |  n  n  n 

18.  When  a  horizontal  beam  exists  only  at  the  ground  surface,  just  one 

set  of  four  eigenvalues  and  two  eigenvectors  have  an  imaginary  part  and  all 

other  sets  (N-l  sets)  do  not  have  such  a  part.  Thus,  equation  25  is  deduced 

to 


(W(x)f  =  [ A (+x) ]  { C . }  +  [A (-x) ]  {C_ }  +  C-  { B (+x)  t  +  C,  { B  C-x)  }  (26 

1  i  l  4 


where,  setting  the  complex  eigenvalue  as  the  first  eigenvalue  (e.g.,  A^  = 


+  iBj  ,  A^  =  ,  A  =  d^  ,  - ),  *B(±x)'  is  expressed  by 


+  d  ^ 

’  B  fix)  }  =  e  *  {cos  fB(x)  {d|j  -  sin  (B^x)  {d^j^ 


(27 


Thus,  the  vertical  displacements  and  shear  forces  induced  in  the  soil  model 
can  be  expressed  in  the  following  forms,  respectively. 


tW(x) }  =  [AB(+x) ] 


Kr ki 


where 


{Ts(x)}  =  [ N] [ AB ' (+x) ]  |c  j  +  [N]  [AB'(-x)]  |c 


[ AB(ix) ]  =  [ [ A(±x) ] I { B (±x) } ] 


[AB '  (±x)  ]  =j^  AB(±x)J 
C1  =  {C^  and  C3  =  {C3} 


Since  the  vertical  displacement  of  the  beam,  W,  (x)  ,  is  located  at  the  first 

b 

value  in  {W(x) }  ,  the  displacement  and  forces  induced  in  the  beam  are 
expressed  by 


[AB(-xl ] 


Wb(x)  -  y  [AB(+x)]|cJ  + 


eb(x)=  5  [AB ’ (+x) ]  | 


;l  *  i;f 


[AB'(-x>] 


M^x)  =  -EI  M  [ AB" (+x) ]  CA  -  El 


[ AB" (-x) ] 


Tb(x)  =  -EI 


[AB"'  (+x)] 


fcll 

C  I 

[M 

-  EI 

I  ~  I 

2! 

lo; 

[ AB" ' (-x) ] 


where 


and  0^  =  vertical  and  rotational  displacements  of  the  beam 
M,  and  T,  *=  moment  and  shear  forces  induced  in  the  beam 


0  =  vector  with  the  size  N  -  1  containing  all  zeros 


[AB*  (±x)  ]  =  AB(±)J 


[AB"(±>]  =  AB ( ± ) 

dx 


[AB,,,(±)]  = 


r  ,3  1 

-At  AB(±x) 
_dx 


Equations  28  and  30  are  combined  in  the  following  manner: 


1C  I  C 

=  [ D (+x)  ]  >  +  [ D (-x)  ] 

c2 )  K: 


T(X)  \  =  [ E(x) 1  1  '  +  [ E ( — x ) ]  J  3 

vx)l  Is)  K) 


wh  ere 


W(x)  =  ' W ( x ) 


T ( x )  =  T  (x) 1  +1 

s  to 


Tb(x) 


Infix)  ] 


( A  B  (  ±  x  )  ] 


[AH'  (ty.)  1 


[ N 1 [AH' (ix) ]  -  El  L |  J  [AB'"  (±x) ] 
[  E  ( i  x )  )  =  - j - 

-El  L  |  [ AB" ( t) ] 
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Responses  due  to  a  line  load 
applied  on  an  infinitely  long  beam 

19.  Vertical  and  moment  loads  are  applied  at  the  location  x  =  0  of  an 
infinitely  long  beam  lying  over  the  ground  surface.  Since  the  responses  must 
diminish  with  distance,  the  expressions  given  in  equations  28,  30,  and  31 
contain  only  the  first  and  second  terms  for  x  ^  0  and  x  ^  0  ,  respectively. 
The  boundary  conditions  at  x  =  0  requires 


s:* 

o 

II 

(w(0)  \ 

\vo)|L 

k(o)i 

(33) 


T(0)  )  _  Jt(0)  )  _  JpI 

vo>k  K(0)|r  w 


where 

Subscript  L  and  R  =  left-hand  and  right-hand  sides  with  respect  to 

the  location  of  the  load  applied 

P  =  load  vector  containing  the  applied  load  at  the 
first  location  of  the  vector  and  zeros  at  other 
locations  of  the  vector 

M  =  applied  moment  load 


Substituting  equation  31  into  equation  33,  the  unknown  constants  are  found  as: 


where  [Y]  =  (N+l)  *  (N+l)  matrix  expressed  by 

[Y]  =  (  (E(+0)  ]  [ D (+0)  ] "" 1  [ D ( — 0)  ]  -  [ E (-0)  ]  ] 
=  [  [ E(+0) ]  -  [ E (-0) ]  (D(-O)]"1  [ D (-0) ]  ] 


(34) 


(35) 


Thus,  substituting  the  constants  expressed  in  equation  34  into  equation  28, 
the  responses  of  the  soil  medium  for  a  line  load  applied  at  x  are  expressed 
by 


20 


(36 


(W(x)  }  =  p(o[ab(-|x  -  c|)]{YjJ  +  m(o[ak(-|x  -  C|)]tYs+,} 

{  Tg  (x)  )  =  P(£)[N][aB'(-|x  -  £()]{Yl}  +  M(0[N][aB'(-)x  -  t  j  >□  ( YN+ £  > 

where  { Y , )  and  =  vector  located  at  the  first  and  N+l  the  columns 

1  N+l 

in  the  matrix  [y]  ,  respectively.  Similarly,  substituting  them  into  equa¬ 
tion  30,  the  responses  of  the  beam  for  a  line  load  applied  at  x  =  £  are 
expressed  by 


V 


U'b(x)  =  P(£)  j„j  [AB  (  -  !  x  -  tD^lYj)  +  M(£)  >  [AB  (- 1  x  - 


£b(x)  =  P(C)  I _}  [aB’(-|x  -  C,')]fY,>  +  M(£)  J  [AB’  (-|x  -  C  ;)]  [y 


N+l  " 


(37 


M^x)  =  -El  [ P(0 


[AB"(-|x  -  C|)]  fy  }  +  M(£) 


[AB’(-Ix  -  £|)]  {yN+1} 


Tb(x)  =  -EI(P(C)  J J  [AB'"  (-jx  -  C  J  >]  ( Y  L  }  +  M(C) I  J  [AB’"  (- |x  -  C 1 )]  { YN+ 1  > 


The  expressions  for  [AB]  ,  [AB"]  ,  and  [AB"’]  are  given  in  Appendix  A. 

Responses  due  to  a  trapezoidal 
load  applied  on  infinitely  long  beam 

20.  A  vertical  trapezoidal  load,  P(£)  =  aC  +  b  ,  is  applied  over  the 

area  x^  ~  x^  of  an  infinitely  long  beam  lying  over  the  ground  surface. 

Integrating  the  responses  due  to  a  line  load  expressed  in  equations  36  and  37 

with  M(E)  =  0  ,  the  following  expressions  for  the  responses  due  to  this  load 

are  obtained: 


21 


The  expressions  for  [D  ]  ~  fD'"  ]  and  [D,  ]  ~  [D'"]  in  equation  38  are  given 

3  3  D  D 

in  Appendix  A. 

Responses  due  to  loads 
applied  on  a  finite  beam 

21.  A  finite  beam  resting  on  a  horizontal  ground  surface  is  now  consid¬ 
ered.  The  soil  medium  is  divided  into  left,  central,  and  right  areas  as  shown 
in  Figure  5.  In  the  central  area,  the  responses  are  obtained  by  superimposing 


Figure  5.  Beam  resting  on 
soil  medium 


L 


LEFT  1  CENTER  ^  RIGHT 
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the  responses  for  an  infinitely  long  beam  and  the  responses  due  to  the  effects 
of  a  finite  length  of  the  beam.  Thus,  the  responses  of  a  system  made  of  a 
finite  beam  and  soil  are  written  by 


{ W ( x )  }c 

{T  (x)} 
s  C 

Wb(x) 

6b(x) 

VX) 

T,  (x) 


{We(x)  }  +  { W1 (x) } 

(Te(x)}  +  { T1 (x) } 
s  s 

wb(x)  +  wb(x) 
0^(x)  + 

M^(x)  +  tt^(x) 
Tf(x)  +  T.1  ( x ) 


I  Soil 


(40) 


Beam 


where  the  superscripts  i  and  e  indicate,  respectively,  the  responses  for 
an  infinitely  long  beam  and  the  responses  due  to  the  effects  of  a  finite 
length  of  the  beam;  and  the  subscript  C  indicates  the  responses  in  the 

central  area.  The  expressions  for  the  responses  of  an  infinitely  long  beam 

are  given  in  equations  37  and  38  for  a  line  load  applied  and  a  trapezoidal 
load  applied,  respectively.  The  expressions  of  the  responses  due  to  a  finite 
length  are  given  in  equations  28  and  30  in  which  the  unknown  constants  are 
determined  for  the  boundary  conditions  at  both  ends  of  the  finite  beam.  The 

expressions  of  the  soil  responses  in  left  and  right  areas  are  given  in 

equation  16. 

22.  The  boundary  conditions  at  both  ends  of  the  central  area  are 


{W(0)}L  =  {W(0))c 


{W(L)}r  =  {W(L))C 


{Ts(0)}l  -  {Ts(0)}c  -  Tb(0)  j  {0} 


{T  (L) }  -  (T  (L) }  -  T  (L) 

s  R  s  b 


M^O)  =  0 


M^L)  =  0 


The  following  relationships  between  {W}  and  {T^}  also  exist: 


{T  (0)  }_  =  [k]  {W(0)L 

S  L  L 


{Ts(L)}r  =  -[k]  {W(L)}r 


where 


[k]  «  [A'(+0)]  (A(+0) ] 


MM 


rj7 


Substituting  equation  42  into  equation  41,  after  substituting  the  first  two 
boundary  conditions  in  equation  41  into  equation  42,  the  boundary  conditions 
given  in  equation  41  can  be  deduced  to 


[k] 

-[kj 


{W(0)  }c 

fW(L)  )c 

V0) 


(T  (0) } 
s  .C 

{T  (L)  }_ 
s  C 

0 


vo) 

Tb(0) 


{0} 

{0} 


(44) 


Mb(L)  =  0 

Thus,  substituting  the  expressions  for  the  displacements  and  forces  into 
equation  44,  the  boundary  conditions  finally  lead  to  the  following 
simultaneous  equations: 


2 (N+l ) 


2 (N+l ) 


<Mn>  : 

i 

[«12i 

(«2I1  i 

1 

[m22] 

(m  iT  i 

3  « 

i 

(m  iT 

4 

{m  >T  ! 

5  ' 

T 

{M  ) 1 
0 

- 

C! 

'h' 

C2 

h 

>  -  - 

C3 

P3 

i  C4  > 

P. 

(45) 


where  { C  }  ,  C  n 
equation  28, 


{S} 


and  C,  are  unknown  constants  defined  in 
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After  determining  the  unknown  constants  from  equation  45  and  substituting  them 
into  equations  28  and  30,  the  responses  due  to  the  effects  of  a  finite  length 
of  the  beam  (e.g.,  {W6}  ,  {T^}  ,  Wb  ’  ®b  *  *^b  ’  an ^  e3uation  40) 

can  be  obtained.  Then,  these  responses  are  added  to  the  responses  of  an 

infinitely  long  beam  (e.g.,  {W*}  ,  { ,  0*  ,  ,  and  T*  in 

equation  40)  to  obtain  the  responses  of  a  finite  beam.  When  a  number  of  loads 

are  applied  simultaneously,  the  resultants  of  those  induced  by  each  of  the 
applied  loads  are  used  for  the  responses  of  an  infinitely  long  beam  which 
appears  in  equations  40  and  46. 

23.  The  above  formulations  are  coded  in  a  computer  program  2PARM. 
Information  of  the  program  is  given  in  Appendix  B. 


PART  IV:  DETERMINATION  OF  THE  PARAMETERS  E'  AND  G' 


Soil  Medium  and  Axial  Normal  Stress 


24.  For  the  idealized  soil  medium  as  shown  in  Figure  1,  the  axial  nor¬ 
mal  stress  in  the  column  and  shear  stress  along  the  faces  of  the  column  are 
assumed  to  be  expressed  by,  respectively, 


(1  - 


o(x,z)  =  F 


v  ) 
s 


s  3w(x,z) 

E  2(1  +  v  )  (1  -  2v  )  3z 
s  s 


and 


(47) 


t(x,z) 


E 

p  - S - 

G  2(1  +  v  ) 
s 


3w(x, z) 
3x 


where 

Eg  and  vg  =  Young's  modulus  and  Poisson's  ratio  of  the  soil  medium 
F  and  F  =  factors  explained  below 

The  above  expressions  without  F  and  F  correspond  to  those  for  a  contin- 

L  Cj 

uous  medium  under  a  completely  restrained  lateral  displacement.  The  factors 
are  included  in  the  expressions  so  that  the  vertical  displacements  of  the 
model  are  close  to  those  of  the  elastic  medium  subjected  to  vertical  surface 
load.  Comparing  equation  1  with  equation  47,  the  parameters  of  the  beam  ele¬ 
ment  in  the  model  are  found  to  be 


(1  -  v  )  E 

p  - ® 5 - 

E  2(1  +  v  )  (1  -  2v  ) 
s  s 

(48) 

E 

p  - 5 - 

G  2(1  -  v  ) 
s 


E' 


25.  The  shape  of  the  surface  displacement  of  the  model  is  controlled  by 

the  value  .  Referring  to  equations  14  and  23,  one  can  find  that  the  value 

X  depends  on  the  ratio  E'/G'  and  thus  F^./F  .  The  surface  displacements 

of  the  model  with  F  /F  =  1  and  2  and  are  shown  in  the  variations  of  Fig- 

r.  ii 

ure  6  for  the  normalized  value.  A  close  agreement  between  the  computed  values 
and  those  obtained  by  the  finite-element  method  (FEM)  can  be  seen  for 
Fg/F  =  2  ,  in  general. 


Determination  of  F  Values 

h 


26.  Under  the  fixed  F  /F  ratio,  the  absolute  displacement  of  the 

h  G 

model  increases  linearly  with  decreasing  F  without  changing  the  displace¬ 
ment  shape  along  the  horizontal  direction.  For  Fp^/F  =  2  ,  the  values  Fp 
are  determined  so  that  the  surface  displacement  of  the  model  is  identical  to 
the  displacement  of  the  continuous  medium  at  the  center  of  a  uniformly  loaded 
area.  Those  values  are  shown  in  Figure  7.  It  is  seen  in  the  figure  that  the 
values  are  rather  independent  of  the  depth  of  the  soil  medium,  but  dependent 
on  a  Poisson's  ratio. 

27.  Using  the  factors  F^  and  F^  above  found  for  a  uniformly  dis¬ 
tributed  load,  the  vertical  surface  displacements  of  the  model  were  computed 
for  a  linearly  varying  distributed  load.  The  values  are  shown  in  Figures  8a, 
b,  and  c,  together  with  those  obtained  by  a  FEM.  Good  agreements  of  the  com¬ 
puted  displacements  with  those  obtained  by  a  FEM  indicates  that  the  factors 

F  and  F  are  rather  independent  of  a  tvpe  of  the  load  applied  on  the 
h  b 

surface . 

28.  From  the  above  studies,  the  approximate  values  of  F  and  F  in 
equation  48  are  found  to  be 

F_  =  0.95  -  0.99  and  F_  =  0.475  ~  0.495  for  v  =  0  -  0.3 
EG  s 

(49a) 


less  than  the  above  values  for  v  =  0.3  -  0.45 


Figure  8a.  Ground  surface  displacement  computed  for  linearly 

varying  load  (v  =  0.0) 


or  given  by 


F_  =  0.96  -  0.49v  +  6.47v2  -  15.76v3 

E  s  s  s 


F 


G 


0. 5F 


E 


(49b) 


The  factors  F  and  F  are  somewhat  independent  of  the  thickness  of  the  soil 
E  Cr 

medium  and  loading  conditions.  Thus,  the  model  parameters  can  be  uniquely  de¬ 
termined  from  equation  48  with  the  factors  defined  in  equation  49,  if  the  elas 
tic  constants  of  the  soil  medium  are  known. 


PART  V:  ASSESSMENT  OF  ACCURACY 


Determining  Factors 

29.  The  accuracy  of  the  computed  values,  using  the  chosen  model,  de¬ 
pends  on  the  way  in  which  the  medium  is  divided  into  a  number  of  horizontal 

layers  and  selection  of  the  values  for  the  factors  F  and  F  .  Figure  9 

L  G 

shows  the  surface  displacements  due  to  a  uniform  load  applied  on  an  elastic 
stratum  which  is  equally  divided  into  various  numbers  of  horizontal  layers  for 
computation.  The  coarser  the  layering,  the  smaller  is  the  displacement  under 
a  loaded  area. 
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Figure  9.  Ground  surface  displacements  for  various 
numbers  of  layering 


30.  A  homogeneous  soil  medium  is  divided  into  a  number  of  horizontal 
layers  as  shown  in  Figure  10.  The  surface  displacements  of  this  medium  were 
computed  for  the  loading  and  are  shown  in  this  figure.  The  load  considered 
here  can  be  produced  by  combining  uniform  and  linearly  varying  loads  that 


LOAD  I 
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a.  Case  considered 


LOAD  II 


b.  Divided  load  and  soil 
medium  for  computation 


Figure  10.  Case  considered  for  assessment  of  present  method 
have  been  previously  used.  The  factors,  F  =  0.95  and  F  =  0.475  ,  were 

h  Cj 

utilized  in  the  analysis  with  the  linear  displacement  approach.  The  displace 
ments  computed  by  this  method  are  very  close  to  those  obtained  by  a  FEM 
(Figure  11). 
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Figure  11. 


Ground  surface  displacements  computed  for  the 
case  shown  in  Figure  10 


31.  The  displacements  of  a  homogeneous  soil  medium  subjected  to  a  sur¬ 
face  load  were  computed  at  various  depths  with  this  working  model.  The 
comparison  of  those  computed  results  with  those  obtained  by  a  FEM  indicate  that 
this  model  can  produce  the  displacement  inside  the  medium  reasonably  well  as 
shown  in  Figure  12. 

Model  Investigation 


32.  The  performance  of  the  present  model  is  investigated  for  the  beams 

bearing  on  a  homogeneous  soil  medium  under  various  conditions.  The  conditions 

considered  are  shown  in  Figure  13.  The  factors  used  are  F  =  0.95  and 

E 

F  =  0.475  and  the  soil  medium  is  divided  into  a  number  of  horizontal  layers 
in  two  different  ways  as  shown  in  Figure  14.  The  responses  of  the  beam  and 
soil  reactions  were  computed  for  vertical  line  loads,  moment  loads,  and  a  uni¬ 
form  load  applied  on  the  beam.  The  computed  results  are  plotted  in  Figures  15 
and  16  for  the  loads  applied  symmetrically  with  respect  to  the  center  of  the 
beam,  in  Figure  17  for  the  loads  applied  symmetrically,  and  in  Figure  18  for  a 
combination  of  various  types  of  the  load.  The  values  obtained  by  a  finite- 
element  method  are  also  plotted  in  the  figures  for  a  comparison.  It  is  seen 
in  those  figures  that  the  present  model  can  accurately  predict  the  beam  re¬ 
sponses  even  with  a  small  number  of  layers. 
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Figure  13.  Case  considered  for  assessment  of 
present  method 
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Figure  14.  Soil  medium  divided  into  number  of  layers 
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Figure  17b.  Responses  computed  for  the  case  considered  in 
Figure  13  (unsymmetric  moment  load,  H  =  100') 


Figure  18.  Responses  computed  for  the  case  considered  in 
Figure  13  (combination  of  mixed  loads,  H  =  100') 


PART  VI:  SUMMARIES  AND  CONCLUSIONS 


33.  This  report  involves  the  development  of  a  two-parameter  model  as  a 
layer  for  the  analysis  of  the  responses  of  a  beam  laying  over  a  ground  sur¬ 
face.  This  mechanical  model  is  made  of  a  number  of  vertical  elastic  columns 
defined  by  the  two-material  parameters  associated  with  axial  force  in  the 
column  and  friction  force  along  the  sides  of  the  column.  The  material  param¬ 
eters  are  defined  by  the  elastic  constants  of  the  medium  (E  and  v  )  and  the 

s  s 

factors  F£  and  F  .  The  factors  are  rather  independent  of  the  thickness  of 
the  median  and  loading  conditions  at  the  surface  of  the  medium,  but  dependent 
on  Poisson’s  ratio.  The  factor  F  is  found  to  be  half  of  the  factor 

u 

=0.5  Fg)  and  the  values  for  F^  can  be  obtained  from  Figure  7  or 
equation  49. 

34.  Using  the  developed  soil  model,  analytical  expressions  for  re¬ 
sponses  of  the  beam,  laying  over  a  ground  surface  and  subjected  to  vertical 
loads,  were  obtained.  These  expressions  are  coded  in  a  computer  program 
2PARM.  Performance  of  the  soil  model  and  the  developed  computer  program  were 
examined  under  various  conditions  by  comparing  the  computed  results  with  those 
obtained  by  a  FEM.  This  study  has  confirmed  that  the  model  and  computer 
program  developed  can  well  predict  the  responses  of  the  beam  bearing  on  the 
surface  of  an  elastic  homogeneous  medium. 
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APPENDIX  A:  "2PARM"  PROGRAM  INFORMATION 
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I.  The  matrices  [AB]  ,  [ AB * ]  ,  [AB"]  ,  and  [AB"']  contain  N+I 

columns  with  a  size  N  in  which  N  is  the  numbers  of  the  layers  in  the  me¬ 
dium.  The  expressions  of  the  N  column  in  these  matrices  are  given  by: 
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first  position  in  the  vectors  given  in  equation  B2. 
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APPENDIX  B:  COMPUTER  PROGRAM  "2PARM1, 


1.  The  computer  program  "2PARM"  is  used  for  the  elastic  response  anal¬ 
ysis  of  a  beam  lying  over  a  ground  surface.  The  formulations  in  the  program 
are  obtained  by  idealizing  the  soil  medium  as  a  new  two-parameter  model  ex¬ 
plained  in  the  text.  The  program  can  handle  the  following  conditions: 

a.  Elastic  soil  defined  by  Young's  modulus  and  Poisson's  ratio. 

b.  With  and  without  elastic  uniform  beam  lying  over  a  horizontal 
ground  surface. 

c.  Horizontally  layered  soil  medium,  including  a  homogeneous 
medium,  underlain  by  a  rigid  base. 

d.  A  number  of  vertically  distributed  and/or  simultaneously 
applied  line  loads  at  any  location  on  the  beam. 

e.  Plane  strain  conditions. 

2.  In  the  analysis  the  soil  is  divided  into  a  number  of  horizontal 
layers  with  soil  properties  linearly  varying  within  a  layer,  and  the  dis¬ 
tributed  load  is  divided  into  a  number  of  trapezoidal  loads. 

3.  The  program  computes  the  following  at  any  location  along  the  ground 
surface,  beam,  and  interface  between  the  layers: 

a.  Shear  and  moment  forces  induced  in  the  beam. 

b.  Vertical  and  rotational  displacements  and  forces  induced  in  the 
beams . 

c.  Vertical  displacements  at  the  interfaces  between  the  soil 
layers. 

d.  Contact  soil  reaction  against  the  beam. 

4.  The  computation  process  is  explained  in  Figure  Bl.  The  following 
are  the  descriptions  of  the  subroutines  in  the  program: 

INPUT  This  subroutine  reads  in  all  necessary  data  that  is  related  to  the 
problem. 

SKMAT  It  forms  the  generalized  stiffness  matrix,  [k]  ,  in  the  equation 
[EI]{WIV}  -  [N] {W"  }  +  [K ] {W }  =  {P}. 

DMAT  It  forms  the  generalized  shear  matrix,  [N]  ,  in  the  equation  above. 

EIGEN  It  computes  the  eigenvalues/eigenvectors  for  soil  regions  where  beam 

is  not  present  by  using  Jacobi  iteration  method.  Ax  ■  x  . 


Bl 


No  Beam 


EIGNC  It  computes  the  eigenvalues/eigenvectors  for  region  where  beam  is 
present.  Ax  =  Bx  . 

OZHES  It  reduces  the  generalized  real  eigen  problem  Ax  =  Bx  to  an  equiv¬ 

alent  problem  with  A  in  upper  Hessenberg  form  and  B  in  upper 
triangular  form,  using  orthogonal  transformations. 

QZIT  It  reduces  A  and  B  as  computed  in  subroutine  OZHES  to  equivalent 

problem  with  A  now  in  quasi-upper  triangular  form  and  B  still  in 
upper  triangular  form. 

OZVAL  A  and  B  from  subroutine  QZIT  is  used  to  compute  quantities  whose 
ratios  give  the  generalized  eigenvalues. 

OZVEC  It  determines  the  eigenvectors  corresponding  to  the  set  of  specified 
eigenvalues  as  found  in  subroutine  OZVAL. 

MATRIX  Formation  of  matrices  using  the  eigenvalues/eigenvectors  as  found  in 
subroutine  EIGEN  and  EIGNC. 

FHMAT  Formation  of  matrices  in  the  solution  expression. 

SDISP  Computes  displacements  for  cases  without  beam. 

SLOAD  This  subroutine  is  called  by  subroutine  SDISP  to  compute  the  solu¬ 
tion  due  to  all  load  cases. 

SPOINT  Computes  solution  due  to  point  loads  on  soil  surface. 

DISTR  Computes  the  solution  due  to  distributed  loads. 

BINFIT  Computes  the  solution  for  point  load  on  infinite  beam. 

BSOLT  Computes  the  unknown  constants  in  solution  expressions  for  the  case 

with  beam. 


BLOAD 


Computes  the  particular  solution  due  to  all  load  cases  on  beam. 


BPOINT 

AINV 

DISPL 

BDISP 


Computes  the  particular  solution  for  point  loads  on  infinite  beam. 

Subroutine  for  real  matrix  inversion. 

Computes  the  solution  for  soil  regions  outside  the  loaded  beam. 

Computes  the  displacement,  moment,  shear,  and  contact  forces  for 
regions  occupied  with  beam. 


Inputs 

5.  The  input  variable  names  are  explained  in  Figure  B2.  The  following 
is  the  input  card  format  for  Program  2PARM: 


Card  Group  I  Identification  Card  -  FORMAT  (20A4) 

1-80  Alphanumeric  identifying  information  to  be  printed  as  titles  on 

the  output 


Card  Group  2  Control  cards  -  FORMAT  (110) 

1-10  N  -  Number  of  soil  layers  (maximum  =  10) 


Card  Group  3 
First  card: 

1  -  10 

Second  card: 


Material  property  information 

Soil  material  property  information  -  FORMAT  (F10.3) 

Poisson's  ratio  of  soil 

Young's  modulus  of  soil  assumed  linearly  varying  with  depth  - 
FORMAT  (2F10.3) 


1-10  El (I)  -  Young's  modulus  at  top  fibre  of  soil  layer  I  (psf) 

11  -  21  E2(I)  -  Young's  modulus  at  bottom  fibre  of  soil  layer  I  (psf) 

Additional  card  must  be  supplied  for  each  soil  layer. 

Third  card:  Soil  layer  thickness  -  FORMAT  (7F10.3) 

1-10  THK(I)  -  Soil  thickness  of  the  Ith  layer  (ft) 

Repeat  for  the  remaining  six  fields  and  continue  on  additional  cards  if 
necessary. 

Fourth  card:  Beam  geometry  and  material  property  -  FORMAT  (3F10.3) 

1-10  BL  -  Length  of  beam  (ft) 
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Figure  B2.  Input  data  notations  used  in  computer  program 

11  -  20  EIBF.AM  -  Young's  modulus  of  beam  (psf) 

21  -  31  ZIBEAM  -  Moment  of  inertia  of  beam  (ft  **  4) 

Card  Group  4  Loading  conditions 

First  card:  Total  number  of  loads  -  FORMAT  (3110) 

I- 10  NPT  -  Number  of  concentrated  point  loads 

II- 20  NMONT  -  Number  of  applied  moments 

21  -  30  NDIRTR  -  Number  of  distributed  loads 

Second  card:  Point  load  magnitude  (downward  positive)  and  their  applied 
locations  -  FORMAT  (2  FI  0.3) 

I  “  10  PT(I)  -  Magnitude  of  point  load  I  (lb) 

II  -  20  DPT(I)  -  Location  of  applied  load  T  (ft) 


Additional  card  must  be  supplied  for  each  point  load. 

Third  card:  Moment  magnitude  (counter  clockwise  is  positive)  and  their 

applied  locations  -  FORMAT  (2F10.3) 

I- 10  CMONT(I)  -  Magnitude  of  applied  moment  I  (ft-lb) 

II- 20  DMT(I)  -  Location  of  applied  moment  I  (ft) 

Additional  card  must  be  supplied  for  each  applied  moment. 

Fourth  card:  Distributed  loads  (assumed  linear  variation;  downward  is 
positive)  and  their  applied  locations  -  FORMAT  (4F10.3) 

I- 10  PI (I)  -  Magnitude  at  near  end  of  distributed  load  I  (lb) 

II- 20  P2(I)  -  Magnitude  at  far  end  of  distributed  load  I  (lb) 

21  -  30  ETA1(I)  -  Near  end  location  of  distributed  load  I  (ft) 

31  -  40  ETA2(I)  -  Far  end  location  of  distributed  load  I  (ft) 

Additional  card  must  be  supplied  for  each  distributed  load. 

Note:  Fields  E  and  F  are  interchangeable.  If  a  number  is  more  suited  to 

one  field,  that  field  may  be  used. 

Card  Group  5  Locations  where  solution  is  desired  -  FORMAT  (3F10.3) 

I- 10  PT1  -  Near  end  point  where  solution  is  to  be  computed  (ft) 

II- 20  PT2  -  Far  end  point  where  solution  is  to  be  computed  (ft) 

21-30  DELTA  -  Regular  interval  between  PT1  and  PT2;  where  the  solution 

is  to  be  computed  (ft) 


Card  Group  6  Modification  factors  -  FORMAT  (2F10.3) 

I- 10  FE  -  Modification  factor  for  axial  forces  in  soil  column 

II- 20  FG  -  Modification  factor  for  shear  forces  in  soil  column 

If  this  card  is  blank,  the  modification  factors  are  computed  from 
Equation  49b. 

6.  Several  problems  may  be  analyzed  in  one  run.  The  data  set  for  each 
problem  commences  with  an  identification  card.  Two  consecutive  blank  cards 
are  required  to  stop  program. 

Outputs 

7.  The  program  prints  the  following  outputs: 

a.  Input  and  generated  data. 

b.  Displacement  (ft),  shear  (k) ,  moment  (k-ft)  and  contact  pres¬ 
sure  (ksf)  at  each  location  as  specified,  card  group  5. 


EXAMPLE  RUN 


INPUT  INFORMATION 


TEST  RUN  1  --  BEAM  WITH  MULTIPLE  LOADING  CONDITIONS 


NUMBER  OF  SOIL  LAYERS,  N 

= 

2 

POISSON'S  RATIO  OF  SOIL 

= 

0.300 

LENGTH  OF  BEAM,  BL 

= 

100.000 

FEET 

YOUNG'S  MODULUS  OF  BEAM 

= 

0. 5160+11 

PSF 

MOMENT  OF  INERTIA  OF  BEAM 

= 

0.100 

FT. **4 

NEAR  END  PT  WHERE  SOLUTION  IS 

DESIRED 

= 

0.000 

FEET 

FAR  END  PT  WHERE  SOLUTION  IS 

DESIRED 

= 

100.000 

FEET 

INTERVALS  BETWEEN  POINTS 

= 

50.000 

FEET 

SOIL  COLUMN  AXIAL  FORCES  MOD. 

FACTOR 

= 

1.000 

SOIL  COLUMN  SHEAR  FORCES  MOD. 

FACTOR 

= 

1.000 

SOIL  LAYER  #  THICKNESS  El ( I )  E2(I) 

1  25.0000  0.4320D+06  0.5760D+06 

2  50.0000  0. 5760D+06  0.5760D+06 


NUMBER  OF  CONCENTRATED  POINT  LOADS  =  2 

POINT  LOAD  (  1)  =  0. 1000D+05  LBS.  AT  X  =  10.000  FEET. 

POINT  LOAD  (  2)  -  0.2000D+05  LBS.  AT  X  =  75.000  FEET. 

NUMBER  OF  APPLIED  MOMENTS  =  1 

MOMENT  (  1)  =  0 . 1000D+06  FT-LBS  AT  X  =  90.000  FEET. 

NUMBER  OF  DISTRIBUTED  LOADS  =  1 

DISTRIBUTED  LOAD  (  1)  : 

PI  =  0 . 5000D+04  LBS.  P2  =  0.1000D+05  LBS. 

ETA1  =  40.000  FEET  ETA2  =  65.000  FEET 


